CU-TP-1153 
FERMILAB-PUB-06-325-T 



The Kaon 5-parameter in 
Mixed Action Chiral Perturbation Theory 

C. Aubin/'Q Jack Laiho,^'!] and Ruth S. Van de Water^fl 

^Department of Physics, Columbia University, New York, NY 10027 
^Theoretical Physics Department, Fermilab, Batavia, IL 60510 
(Dated: February 1, 2008) 

Abstract 

We calculate the kaon B-parameter, Bk, in chiral perturbation theory for a partially quenched, 
mixed action theory with Ginsparg- Wilson valence quarks and staggered sea quarks. We find that 
the resulting expression is similar to that in the continuum, and in fact has only two additional 
unknown parameters. At one- loop order, taste-symmetry violations in the staggered sea sector 
only contribute to flavor-disconnected diagrams by generating an 0{a^) shift to the masses of 
taste-singlet sea-sea mesons. Lattice discretization errors also give rise to an analytic term which 
shifts the tree-level value of Bk by an amount of 0{a^). This term, however, is not strictly due to 
taste-breaking, and is therefore also present in the expression for Bk for pure G-W lattice fermions. 
We also present a numerical study of the mixed Bk expression in order to demonstrate that both 
discretization errors and finite volume effects are small and under control on the MILC improved 
staggered lattices. 
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I. INTRODUCTION 



Lattice quantum chromodynamics allows nonperturbative calculations of low-energy 
QCD quantities from first principles. Until recently, limited computing resources and the 
inability to include quark loop effects have prevented lattice calculations from achieving 
realistic results. In the past few years, however, lattice simulations of both light meson and 
heavy-light meson quantities with dynamical staggered quarks have shown excellent numeri- 
cal agreement with experimental results |l| . These have included both post-dictions, such as 
;he pion and kaon decay constants [2^, and predictions, as in the case of the Be meson mass 
Such successes demonstrate that many of the systematic uncertainties associated with 
lattice simulations are under control, and therefore give confidence that lattice simulations 
can reliably calculate quantities that cannot be accessed experimentally. One of the simplest 
quantities of phenomenological importance that can only be calculated using lattice QCD 
is the kaon B-parameter, B^. Thus a precise measurement of B^ is an important goal for 
the lattice community.^ 

Bk parameterizes the hadronic contribution to K^—K^ mixing; it therefore plays a crucial 
role in extracting information about the CKM matrix using measurements of the neutral 
kaon system. In particular, the size of indirect CP violation in the neutral kaon system, 
when combined with a numerical value for Bx-, places an important constraint on the 
apex of the CKM unitarity triangle [?, 8]. Because ex is well known experimentally [9|], the 
dominant source of error in this procedure is the uncertainty in the lattice determination 
of Bk- It is likely that new physics would give rise to CP- violating phases in addition to 
that of the CKM matrix; such phases would manifest themselves as apparent inconsistencies 
among different measurements of quantities which should be identical within the standard 
CKM picture. Thus a precise determination of B^ will help to constrain physics beyond 
the standard model. 

Because lattice simulations with staggered fermions can at present reach significantly 



lighter quark masses than those with other standard discretizations [10|, lUl], calculations 
of weak matrix elements using the available 2+1-flavor Asqtad staggered lattices appear 
promising. Unfortunately, however, one pays a significant price for the computational speed 



^ Promising calculations of Bk including dynamical quark effects are currently in progress using both 
improved staggered fermions 4, 5] and domain-wall fermions 
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of staggered simulations: each flavor of staggered quark comes in four species, or "tastes." 
In the continuum limit, these species become degenerate and can be removed by taking the 
fourth root of the fermion determinant. In practice, however, one must take the fourth root 
of the quark determinant during the lattice simulation in order to remove the extra tastes; 
thus it is an open theoretical question whether or not one recovers QCD after taking the 
continuum limit of fourth-rooted lattice simulations. Recently several papers have appeared 



13|, 
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16| . Although the issue 



addressing the validity of the "fourth-root trick" |12l . 
has not yet been resolved, there are indications that the fourth root does not introduce 
pathologies when taking the continuum limit of the lattice theory. We recommend a recent 
review by Sharpe, Ref. [l3], as a clear summary of the current status of the fourth-root 
trick. In this paper we assume the validity of the fourth-root trick. Even working under 
this assumption, however, the additional tastes introduce complications to staggered lattice 
simulations. The degeneracy among the four tastes is broken by the nonzero lattice spacing, 
a, and results in discretization errors of (9(a^). Thus one must use functional forms calculated 
in staggered chiral perturbation theory (SXPT), in which taste-violating effects are explicit. 
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2l|. 



to correctly extrapolate staggered lattice data [18 . 

Staggered chiral perturbation theory has been used to successfully fit quantities such as 
/,r and /d, for which the SXPT expressions are simple and taste-symmetry breaking pri- 

31. In the 



marily enters through additive corrections to the pion masses inside loops |2C 
case of Bki however, the additional tastes also make the matching procedure between the 
lattice AS* = 2 effective four- fermion operator and the desired continuum operator more 
difficult. The latticized version of the continuum Bk operator mixes with all other lattice 
operators that are in the same representation of the staggered lattice symmetry group [23|] 
- including those with different tastes than the valence mesons. Current staggered calcula- 
tions only account for mixing with operators of the correct taste and only to 1-loop order in 
as', this results in a 20% uncertainty in Bk ^. In order to achieve a precise determination 
of Bk with staggered fermions, one must either perform the lattice-to-continuum matching 



nonperturbatively using a method such as that of the Rome- Southampton group 24 1 or one 
must include the effects of the extra operators in the SXPT expression for Bk 25|. Thus 
far the large number of staggered operators has prevented matching calculations beyond 
1-loop order H], although nonperturbative matching including all relevant staggered 
operators is, in principle, possible. The effects of this truncated lattice-to-continuum opera- 
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tor matching can be included in an extended version of SXPT, but the resulting expression 
for Bk has many undetermined fit parameters, only a few of which are already known from 
measurements of other quantities. Therefore use of this expression may prove just as difficult 
as implementing fully nonperturbative matching. 

The calculation of weak matrix elements such as with Ginsparg- Wilson (G-W) quarks 



28l | , on the other hand, is theoretically much cleaner than that with staggered quarks. This 



is because in the massless limit, G-W quarks possess an exact chiral symmetry on the 
lattice and do not occur in multiple species. Although in practice, lattice simulations use 
approximate G-W fermions, the degree to which chiral symmetry is broken in simulations 
can be controlled either by the length of the fifth dimension in the case of domain-wall 
quarks 3 '^^ through the degree to which the overlap operator is realized in the case 



of overlap quarks 



32 



33|. Consequently, while the AS* = 2 lattice operator still mixes 



with wrong-chirality operators, there are significantly fewer operators than in the staggered 
case, and nonperturbative renormalization can be used in the determination of Bk- As an 
additional benefit, the approximate chiral symmetry also ensures that the appropriate chiral 
perturbation theory expression for use in the extrapolation of Bk lattice data is continuum- 
like at next-to- leading order (NLO). Lattice simulations with G-W fermions, however, are 
10 to 100 times more computationally expensive than those with staggered fermions with 
comparable masses and lattice spacing, and thus are unfortunately not yet practical for 
realizing light dynamical quark masses. 

A computationally affordable compromise is therefore to calculate correlation functions 
with Ginsparg- Wilson valence quarks on a background of dynamical staggered gauge con- 
figurations. This "mixed action" approach combines the advantages of staggered and G-W 
fermions while not suffering from their major disadvantages, and is well-suited to the cal- 
culation of the weak matrix element Bk- By using staggered sea quarks and G-W valence 
quarks one can better approach the chiral regime in the sea sector while minimizing operator 
mixing and allowing the use of nonperturbative renormalization. Additionally, because the 
MILC staggered lattices with 2-1-1 flavors of dynamical quarks are publicly available and 
offer a number of quark masses and lattice spacings [ll|, one can perform unquenched 
three-flavor simulations at the same cost as quenched G-W simulations. Mixed action simu- 
ations have already been successfully used to study quantities of interest to nuclear physics 
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351 ]. thus we expect that a similar method can be used to calculate Bk- 
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One might be concerned that the use of a mixed action could introduce new theoretical 
complications into the lattice determination of Bk- Although the mixed staggered sea, G-W 
valence lattice theory reduces to QCD in the continuum limit, at nonzero lattice spacing, 
it is manifestly unphysical in that it violates unitarity. Thus, in order to extract physical 
QCD quantities from mixed action simulations, it is essential that one can correctly describe 
and remove contributions due to unphysical mixed action effects from quantities such as Bk 
using the appropriate lattice chiral perturbation theory. It has not been rigorously proven 
that the mixed action chiral perturbation theory developed in jsol is the correct chiral 
effective theory for mixed G-W, staggered simulations in which the fourth-root of the quark 
determinant is taken in the sea sector. Nevertheless, Ref. |l2] showed that, given a few 
plausible assumptions, staggered chiral perturbation theory is the correct chiral effective 
theory for describing the pseudo-Goldstone boson sector of rooted staggered simulations. A 
similar line of reasoning should hold for mixed action chiral perturbation theory. Assuming 
so, mixed action lattice simulations can be used to correctly calculate quantities involving 
pseudo-Goldstone bosons (such as Bk) in QCD. More comphcated quantities, however, 
should be considered on a case-by-case basis. 

In this paper we calculate Bk in XPT for a lattice theory with G-W valence quarks and 
staggered sea quarks. We present results for a "1+1+1" theory in which m„ 7^ 7^ 
in the sea sector, for a "2+1" theory in which m„ = 7^ in the sea sector, and for 
a "full QCD" -like expression in which we tune the valence- valence meson masses equal to 
the taste-singlet sea-sea meson masses. (We emphasize, however, that these expressions 
only truly reduce to QCD when mgea = "^valence and the lattice spacing a 0.) These 
expressions will be necessary for the correct chiral and continuum extrapolation of mixed- 
action Bk lattice data. We find that the expression for Bk in mixed action XPT has only 
two more parameters than in the continuum. The first coefficient multiplies an analytic 
term which shifts the tree-level value of Bk by an amount of 0{a^)] this term is also present 
in the case of pure G-W lattice fermions. The second new parameter shifts the mass of 
the taste-singlet sea-sea meson (which only appears inside loop diagrams) by an amount of 
0(0^). This mass-splitting has already been separately determined, however, in the MILC 
spectrum calculations so we do not consider it to be an unknown fit parameter. Therefore, in 
practice, the chiral and continuum extrapolation of mixed action Bk lattice data should be 
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no more complicated than that of domain-wall lattice data. A numerical analysis using the 
taste-breaking parameters measured on the MILC coarse lattices (a ^ 0.125 fm) shows that 
the size of non-analytic discretization errors should be less than a percent of the continuum 
value of Bk over the relevant extrapolation range. In addition, we find that finite volume 
effects in are also small, and are of 0{1%) for the lightest pion mass on the MILC coarse 
ensemble. 

This paper is organized as follows. We review mixed action chiral perturbation theory 
(MAXPT) in Sec. M In Sec. [ml we calculate Bk to next-to-leading order in MAXPT. This 
is divided into four subsections: we first present the spurion analysis necessary to map the 
quark-level B^ operator onto an operator in the chiral effective theory in Sec. IIII At calculate 
the 1-loop contributions to B^ at NLO in Sec. IIIIB[ determine the analytic contributions 
to Bk at NLO in Sec. IIII CI and finally present the complete expressions for Bk at NLO 
in MAXPT in Sec. IIIIDI Next, in Sec. IIV| we estimate the numerical size of both taste- 
symmetry breaking contributions and finite volume effects on the existing MILC ensembles 
using the resulting mixed action XPT formulae. Finally, we conclude in Sec. |Vl 



II. MIXED ACTION CHIRAL PERTURBATION THEORY 



In this section we review the leading-order mixed action chiral Lagrangian, first deter- 



mined in Ref. [SGj, and discuss some of its physical consequences for the pseudo-Goldstone 
boson sector. 

We consider a partially quenched theory with A^ai Ginsparg- Wilson valence quarks and 
Asea staggered sea quarks. Each staggered sea quark comes in four tastes, and each G-W 
valence quark has a corresponding bosonic ghost partner to cancel its contribution to loop 
diagrams. For example, in the case Aval = 2 and A'sea = 3 , the quark mass matrix is given 
by 

M = diag( m^, m„, m^, m^^ rrid, mq, rris, rris, rUs, nis -, m^^m^ m^JUy) , ( 1 ) 

sea valence ghost 

where we label the dynamical quarks by u, d, and s and the valence quarks by x and y. Near 
the chiral and continuum limits, the mixed action theory has an approximate SU{4Nsea. + 
Avail A^vaOi ® 5'{7(4Asea + A^vai| Avai)_R graded chiral symmetry. In analogy with QCD, we 
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assume that chiral symmetry spontaneously breaks to its vector subgroup, 



5t/(47V,ea + A^val|iVval)L ® SU{4N,,, + N,^\N,^)r ^ SU{4N,,, + iVval|A^val)y, (2) 

and gives rise to (4iVsea + 2A^vai)^ — 1 pseudo-Goldstone bosons (PGBs). These PGBs can 
be packaged in an SU{4:Nsea. + -^vail-^vai) matrix: 

f u 7r+ x+ \ 



exp 



2i$ 



dy 



TT- D Qrf, Q 

K- S Qs. Qsy 



Qux Q\x Qlx ^jjx 
Qiy Qdy Qly ^ ^ ^ly ^\y 



(3) 



where / is normahzed such that /^r 
PGBs, each of which comes in sixteen tastes. For example, 

16 



V ^yx Ryy P ^ / 

131 MeV. The upper-left block of $ contains sea-sea 



(4) 



6=1 



where the Euclidean gamma matrices 



(5) 



are the generators of the continuum SU{A) taste symmetry (^7 is the 4x4 identity matrix). 
The fields in the central block are the flavor-charged (P"^ and P") and flavor-neutral {X and 
Y) valence-valence PGBs, while those in the lower-right block with tildes are the analogous 
PGBs composed of only ghost quarks. Finally, the off-diagonal blocks contain "mixed" 
PGBs: those labelled by i?'s arc composed of one valence and one ghost quark, while those 
labelled by Q's are composed of one valence and one sea quark. We do not show the mixed 
ghost-sea PGBs exphcitly; their locations are indicated by eUipses. 
Under chiral symmetry transformations, E transforms as 



(6) 



The standard mixed action chiral perturbation theory power-counting scheme is 

Ppgb/K ~ ^?/^QCD ~ o^Aqcd , (7) 
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so the lowest-order, 0{ppQ^,mq, a^), mixed action chiral Lagrangian is 

L = ^Str(9^S9^St) - /^Str(SMt + MSt) + (y^ + % + \]y) , (8) 
o 4 

where Str indicates a graded supertrace over both flavor and taste indices and fi is an 
undetermined dimensionful parameter of (9(Aqcd)- The leading-order expression for the 
mass-squared of a valence-valence PGB is identical to that of the continuum because the 
chiral symmetry of the valence sector prevents an additive shift due to lattice spacing effects: 

mly = ^{m^ + my). (9) 

Us and Wg comprise the well-known staggered potential and come from taste-symmetry 
breaking in the sea quark sector [3]. Us splits the tree-level masses of the sea-sea PGBs 
into degenerate groups: 



ff',t 



fi{mf + nif/) + a^At, (10) 



where Aj is different for each of the S'0(4)-taste irreps: P, V, A, T, I. In particular, Ap = 
because the taste-pseudoscalar sea-sea PGB is a true lattice Goldstone boson. The mixed 
action Lagrangian contains only one new operator, and thus one new low-energy constant, 
as compared to the staggered chiral Lagrangian: 

Uy = -CMixStr(r3Sr3St), (11) 

where 

T"3 = Psea " Pval = diag(/sca ® /taste, -/val, -/val)- (12) 

This operator links the valence and sea sectors and generates a shift in the mass-squared of 



a mixed valence-sea PGB: 



16C, 



m 



j-^ ^/i(m/-Fm^) + a^AMix, Awix = • (13) 

Although the parameter Anix has not yet been calculated in mixed action lattice simulations, 
it can, in principle, be determined by calculating the mass of a mixed valence-sea meson on 
the lattice. As in any partially quenched theory, the mixed action theory contains flavor- 
neutral quark-disconnected hairpin propagators which have double pole contributions. The 
only flavor-neutral propagators that appear in the expression for are those with two 
valence quarks; these have the following form: 

Gxyiq) = -^^ + DMq\ (14) 
q -\- nix 



where 

is the disconnected (hairpin) contribution. Note that the sea-sea PGBs in the above ex- 
pression are all taste singlets because the valence quarks do not transform under the taste 
symmetry. 



III. Bk at nlo in mixed action chiral perturbation theory 

In this section we outline the calculation of Bk in mixed action XPT. We divide it into 
four subsections. We first determine the operators that contribute to Bx in the mixed 
action chiral effective theory using a spurion analysis in Sec. IIII A[ In Sec. IIII Bl we outline 
the 1-loop calculation of B^-, and in Sec. IIII CI we follow this up with an enumeration of the 
corresponding analytic terms. Finally, the complete NLO results are presented in Sec. IIIIDI 



A. Bk Spurion Analysis 

The spurion analysis for Bk in the mixed action case is similar to that in the continuum 



371 ]: we will point out differences when they occur. 



In continuum QCD, Bk is defined as a ratio of matrix elements: 



Bk = j—. (16) 



' vac 



The numerator in the above expression measures the hadronic contribution to neutral kaon 
mixing: 

Mk = (k'IOkIK'), (17) 
Ok = [s-f,{l--f,)d][s^,{l--f5)d], (18) 

where we have dropped color indices in Ok because both color contractions give rise to the 
same operators in the chiral effective theory. Ok is an electroweak operator which transforms 
as a {27l, 1_r) under the standard continuum chiral symmetry group. The denominator in 
Eq. f|T6l) is the same matrix element as in the numerator evaluated in the vacuum saturation 
approximation: 



MT = ^{K \s^,{l - 75)c^|0)(0|s7^(l - ^,)d\K'), (19) 



so that Bk is normalized to be of 0{1). 

In the mixed action theory, we define Bk in an analogous manner, except that both the 
external kaons and the operator Ok must now be composed of valence quarks: 



Of = [y7,{l - l,)x][y^,il - 75)x]. (20) 
We can rewrite this operator as follows: 

Of = ^qdl, ® PyMqLil^. ® PyM , (21) 

where the subscript "L" on the quark field indicates the left-handed projection and the 
matrix By^. is a 4:Nsea + 2A^vai matrix in fiavor space that projects out the yx component 
of each quark bilinear {[Byx]ij = 5j,i45j,i3 for Nsca = 3, A^vai = 2). If we promote By^ to 
a spurion field, Bk, which can transform under the mixed action chiral symmetry group 
5f/(4A/'sea + A^vailA^vaOi ® SU{AN,^^ + A^vai | A^vai)i?, then Bk must transform as 



LBkL^ , (22) 



so that Eq. fl2T|) remains invariant. One cannot build a chirally invariant operator out of 
E and the spurion field Bk without derivatives, but one can build two such operators at 

OiplcB)-- 

J2 Str[Sa^StFK]Str[Sa^StF;,], (23) 

J2stT[w,j:^BKJ:d,j:^BK]. (24) 

It turns out however, that these two operators are equivalent when one demotes the spurion 
Bk to the matrix By^. Thus we are left with a single chiral operator: 

O^K = J2 Str[S9^StP5,]Str[Sa^Stp^J . (25) 



37|, but S contains more 



This operator is identical in form to the continuum Bk operator 
fields and the standard trace has been promoted to a supertrace. 

Because this operator is of ^(ppgb)) operators of 0{mq) and C(ci^), if present, could 
also potentially contribute at the same order in XPT. Recall that the quark mass matrix, 
when promoted to a spurion field, must transform as M ^ LMi?^; thus we cannot form a 
chiral operator with a single power of M and two powers of Bk that is invariant under the 
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chiral symmetry. This is, of course, to be expected because the chiral symmetry of the G-W 
valence sector and the U{1)a symmetry of the staggered sea sector are sufficient to prevent 
any new operators involving the quark mass matrix at leading order in the chiral expansion. 
New operators of O(a^) that are not present in the continuum can also potentially appear 
and contribute to B^. As discussed in Ref. [25|, such operators arise in two distinct ways: 
mixing with higher-dimension operators and insertions of four-fermion operators from the 
action. We demonstrate that these operators do not introduce taste-symmetry breaking, 
and therefore give rise to the same kinds of analytic terms as in the pure G-W case. 

Although the Bk lattice operator is dimension 6, at the level of the Symanzik effective 
theory, it maps onto all continuum effective operators of dimensions 6 and higher that re- 
spect the same lattice symmetries. Operators of dimension 7 and 8 are explicitly suppressed 
relative to the dimension 6 Bk Symanzik effective operator by powers of a and a^, respec- 
tively, and can therefore be mapped onto chiral operators that may contribute to Bk at 
NLO. Because, however, the lattice symmetry group includes taste transformations under 
which the valence quarks are singlets, only dimension 7 and 8 operators composed of four 
valence quarks can possibly respect the same lattice symmetries as the Bk lattice operator. 
Moreover, the chiral symmetry of the G-W valence quarks in the Bk operator prohibits 
strictly valence dimension 7 four-fermion operators. Thus we need only consider dimension 
8 Symanzik effective operators for the Bk operator which contain four valence quarks. For- 
tunately we need not enumerate all possible dimension 8 quark-level operators of O^a^p^) 
and 0{a?"w?) in order to determine all possible chiral operators of O^a^p^Q^) and 0{a?mq) 
onto which they map. Because the Bk lattice operator has a L-L chiral structure, and 
chiral symmetry is respected by the valence sector of the lattice theory, it can only mix with 
higher-dimension operators that also have a L-L structure. Consequently, these dimension 8 
four-fermion operators will generate the same spurious as the dimension 6 Bk operator, and 
thus lead to the same chiral operator as in Eq. ( l25l) . but with an additional undetermined 
coefficient of (9(a^). In general, this new coefficient just produces an unknown shift of O{o^) 
to the original 0{1) coefficient of the O^. Thus it will not lead to any new functional forms 
in the expression for Bk in MA^PT, only additional contributions from the leading-order 
Bk operator that of are higher-order in the MA^PT power-counting. In particular, at NLO 
in mixed action xPT, it will simply lead to an 0{a'^) correction to the tree-level value of Bk-, 
which we can absorb into an analytic term. We emphasize that, although this contribution 
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is proportional to a^, it is not due to taste-symmetry breaking in the staggered sea sector. 
Because it arises from strictly valence four-fermion operators, it is also present in simulations 
with pure G-W lattice fermions. 

New operators that contribute to at O(a^) can also be produced by inserting a 
dimension 6 O(a^) operator from the Symanzik action into the four-fermion operator. 
A method for combining four-fermion operators at the chiral level was developed in Ref. 21 1 



for the purpose of determining the NLO staggered chiral Lagrangian. This method was 
later used in Ref. 25(] to enumerate the operators that arise from insertions of the staggered 
action that contribute to Bk with both staggered sea and valence quarks. In the case of 
the mixed action theory, we must consider insertions of four-fermion operators with only sea 
quarks, four-fermion operators with only valence quarks, and four-fermion operators with 
both sea and valence quarks. 

Let us first consider insertions of operators with only sea quarks, as all of the work has 



essentially been done in Ref 251]. Staggered four-fermion operators can be made invariant 



under arbitrary chiral symmetry transformations by introducing six pairs of spurion fields: 

Fl®Fl^ LFlL^ ® LFlL^ , Fr^Fr^ RFrR^ ® RFrR^ , 
Fl®Fr^ LFlL^ ® RFrR\ Fl®Fl-^ LFlR^ ® LFlR\ 
Fr®Fr-^RFrU®RFrL\ Fl^Fr^ LFlR' ®RFrL\ (26) 

where the two separate spurious in each pair correspond to the two separate quark bilinears 
in each four-fermion operator. The taste-breaking spurious in Eq. fl26l) can be combined 
with the Bk spurion F^-, which transforms as in Eq. (122|) . to produce all of the generic 
chiral structures in Table [H^ In order to turn these structures into operators that contribute 
to Bki one must ultimately replace the spurion fields with constant values. In the case of 
the staggered theory, because the constituent staggered bilinears may have nontrivial tastes, 
these spurious become taste matrices, S,i ® ^i, which are diagonal in flavor space. In the 
case of the mixed action theory, because only the staggered sea quarks carry taste quantum 
numbers, these spurious become taste matrices multiplied by projectors onto the sea sector, 
i.e. — > C.iVsea.^C.i'Psea.- For example, consider a particular operator in the mixed action 



Note this table is simply the first column of Table III in Ref. 25 1 
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Generic Chiral Structure 
Str{FK^FR^^FK^F'j^^^) + p.c. 
Str(F^SF^jSt)Str(FKSFKSt) +p.c. 
Str{FKFLE^FK^FR)+p.c. 
Str(F^i^St)Str(Fi^Si^) +p.c. 
Str{FKFL^^FKFLJ:^)+p.c. 
Str(FA'FLSt)Str(Fi^FLSt) + p.c. 

TABLE I: Mesonic operators corresponding to insertions of four-fermion operators from the 
Symanzik effective action; tliese generic structures apply to operators with only sea quarks, to 
those with only valence quarks, and to mixed operators with both sea and valence quarks. Fx 
comes from the Bk quark-level operator and will ultimately be set equal to the projector Pyx- The 
remaining spurious come from four-fermion operators, and must be set equal to different matrices 
depending upon the four-fermion operator under consideration; the specific details are discussed 
in the text. The notation "p.c." indicates the parity-conjugate of the previous operator. 

theory coming from the first spurion combination in Table [H 

StriFnEFKl^^FREFKJ:^) + p.c. Str{C,rse.^Pycc^^^5Vse.^Py.^^) + p.c. (27) 

Because this operator is already of (9(a^), it can only contribute to Bx at NLO through 
tree-level diagrams; thus it can only contain two pion fields, which are insufficient to separate 
all of the projectors onto the sea sector from all of the projectors on the valence sector. In 
fact, it is easy to show that none of the chiral operators arising from insertions of staggered 
four-fermion operators actually contribute to at NLO (although they may at higher 
orders) because their contributions vanish identically due to the fact that VseaPyx = 0. Thus 
we do not show all of their expressions here. 

We next consider insertions of mixed four-fermion operators in which one quark bilinear 
contains staggered sea quarks and the other contains GW valence quarks. Such mixed 
operators can be made invariant under arbitrary chiral symmetry transformations with the 



following four pairs of spurion fields [36|: 



Fl®Fl^ LFlL^ ® LFlL^ , Fr^Fr^ RFrR^ ® RFrR^ , 
Fl^Fr^ LFlL^ ® RFrR\ Fr®Fl-^ RFrR^ ® LFlL\ (28) 
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where, in this case, we have distinguished between the first spurion in the pair (which 
corresponds to the sea bihnear) and the second spurion (which corresponds to the valence 
bihnear). This is necessary because the sea spurion will ultimately be replaced with a 
projector onto the sea sector, while the valence spurion will become a projector on the 
valence sector. By comparing Eqs. fl28p and (12^ . one can see that the mixed spurion fields 
are actually a subset of the staggered spurion fields, although they must be replaced with 
different matrices in order to produce operators that contribute to Bk- Consequently, the 
only possible combinations of the spurion with the mixed spurious are already given 
in the upper panel of Table [H Although the mixed case is clearly similar to the staggered 
one, let us nevertheless consider the example of the operator arising from the first generic 
structure in Table [H 

StriFREFK^^FnEFKT.^) + p.c. "'^'^ StT{Vse.^Pyx^^V,.i^Py,^^) + p.c. (29) 

This operator does not contribute to B^ at NLO for the same reasons as the previous 
example, and neither does the other operator arising from insertions of mixed valence-sea 
four-fermion operators. 

Finally, we consider insertions of purely valence four-fermion operators. These operators 
lead to the same spurion fields as the mixed operators [36], and thus to the same chiral 
forms in the upper panel of Table [B The only difference is that both spurion fields must 
be replaced with projectors onto the valence sector. Consequently, as in the case of the 
previous example, the new chiral operators do not contribute to Bk at NLO. In summary, 
although many MA^PT operators of 0{a^) arise from insertions of four-fermion operators 
in the Symanzik effective action, none of them contribute to Bk at NLO for the mixed G-W, 
staggered lattice theory. 

The previous analysis holds for G-W valence quarks, which have perfect chiral symmetry. 
On the lattice however, G-W quarks are often approximated as domain-wall quarks, which 
have a small amount of residual chiral symmetry breaking due to the finite size of the fifth 
dimension. This chiral symmetry breaking is parameterized by the residual mass, mres, 
which is a measure of how far the left- and right-handed components of the quarks extend 
into the fifth dimension. These effects can be readily added to the chiral theory, as seen in 



Refs. m 



391], by adding the following mass-like term to the chiral Lagrangian: 



^^DWF ^ _/fZ!str (Sl]t + , (30) 
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where Q is a spurion which transforms as the mass matrix transforms, and in the end we 
set Q = rures X I- This leads to the famihar expression for the tree- level mass of a PGB 
composed of two domain-wall quarks: 

mly^li{mx + my + 2m^es)- (31) 

Clearly this term will not contribute at leading order to B^, since the Q spurion transforms 
in the same manner as the mass spurion, and the mass term did not contribute at this order. 
Consequently, one may simply shift the valence quark masses by mvai — irivai + rrires in the 
results of this paper to transform them into expressions that apply to lattice simulations 
with domain-wall valence quarks and staggered sea quarks.^ 

Finally we must consider the fact that, while one constructs a lattice operator to corre- 
spond to a continuum operator with a particular spin, once on the lattice, this operator is 
allowed to mix due to gluon exchange with operators that correspond to other continuum 
spin structures. This comes from the fact that the symmetry group on the lattice is not the 
5*0(4) group of Euclidean rotations: it has been broken down to the subgroup of hypercubic 
rotations. Lattice operator mixing patterns can become especially complicated when one 
introduces staggered quarks because now the desired lattice operator can not only mix with 
other operators with incorrect spins, but also those with incorrect tastes. Fortunately, in 
the mixed action theory that we consider here, the symmetry of the G-W valence sector 
is sufficient to prevent mixing between the lattice Bk operator and new operators with 
nontrivial taste structure. Because the valence quarks in the four-fermion operator do 
not transform under taste symmetry, the Bk operator is clearly not in the same lattice 
symmetry irrep as any taste-violating four-fermion operators. Furthermore, in the case of 
pure G-W valence quarks, the B^ operator cannot mix with operators of the wrong chiral- 
ity. In realistic simulations with domain- wall valence quarks, however, the desired lattice 
Bk operator with spin structure VV + AA mixes with four other operators which do not 
have the same VV + AA spin structure (TT, VV - AA, SS + PP, and SS - PP). This 

^ In applying xPT expressions to lattice simulations with domain-wall valence quarks and staggered sea 
quarks, it is also important to remember that the domain- wall valence and staggered sea quark masses 

are renormalized differently. Consequently, if one wishes to use the bare domain-wall lattice Dirac mass 
parameter and the bare AsqTad staggered lattice mass in mixed action xPT expressions, one must allow 
the parameter /x which relates the quark masses to the pion mass-squared to be different in the valence 
and sea sectors. 
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contamination though is suppressed by two factors of the residual mass [40j, and the effect 
is small, so long as the residual mass is small. Although the effect may be small, it could be 
non- negligible, but then it can be removed nonperturbatively using the standard method of 



Rome-Southampton 



24|. 



To summarize, the result of this spurion analysis is that the leading order operator that 
contributes to Bk in the mixed action case is simply the continuum operator naively gener- 



alized to the mixed action theory. This is simpler than the full staggered case [25|] , in which 
many new operators appeared at leading order due to taste-symmetry breaking. Neverthe- 
less, taste-symmetry breaking will still enter the mixed action calculation of Bk through 
the masses of PGBs in loop diagrams. In addition, taste-breaking operators will generate 
operators of next-to-leading order (NLO) in the chiral effective theory that will contribute 
to analytic terms; we will discuss these in Sec. IIII CI 



B. Contribution of at 1-Loop 

Recall from Eq flTB]) that the kaon B-parameter is defined as the ratio Aixf-M^K'^- At 
tree-level, 

(^)'°-°- 

Because all higher-order contributions to Bk are identically zero in the limit of massless 
quarks, this expression defines the B-parameter in the chiral and continuum limits. 

At 1-loop, the K'^ — 'K^ matrix element receives contributions from the diagrams shown 
in Fig. [H where we have specified the location of each of the two left-handed currents in the 
chiral operator O^.^ This factorization of the operator is useful because the calculation of 
the kaon matrix element can then be separated into two pieces, the 1-loop corrections to Jk 
and the 1-loop corrections to Bk in which we are interested. In terms of the contributions 
from the diagrams in Figure [U the Bk matrix element can be written as 

Mk = lBofm%{l + X[Figs. l(b)-(c)]} + X'[Figs. l(d)-(f)]}, (33) 

Note that this factorization is only possible at leading order. At higher orders, receives contributions 
from operators which are not simply products of left-handed currents. 
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(a) (b) 



^ XI iL 

(d) (e) (f) 

FIG. 1: Tree-level and 1-loop contributions to Mk- The circle represents a vertex from the LO 
staggered chiral Lagrangian. Each square represents an insertion of one of the two left-handed 
currents in and "changes" the quark flavor from d <-> s. 

where X and X' indicate the results of specific diagrams and m'^y is the 1-loop kaon mass 
squared. At 1-loop, the form of TW^'^ is simple: 

MT = lmlyf%, (34) 

where m'^y and fxy are the 1-loop corrected values. It is clear that diagrams (b) and (c) in 
Fig. [Ufactorize - the left-half of each diagram is the 1-loop correction to fx, while the right 
half is just / at tree-level. Mathematically, 

X[Figs. l(b)-(c)] = 2^, (35) 

where the factor of two comes from the fact that the loop can appear on either leg. Diagrams 
Wih) and[T](c) therefore turn the leading order into the 1-loop f^y in Eq. fl33|l : 

Mk = iBoflmly + X'[Figs. l(d)-(f)], (36) 

such that Bk at one loop only depends on diagrams [T](d)-(f): 

„.-i«P „ , 3 Jf'[Figs. l(d)-(f)] 

Figure [2] shows the quark flow diagrams that correspond to the meson diagrams in 
Figs. [T](d)-(f). It is interesting to note that the only place where sea quarks appear in 
these diagrams is in the disconnected hairpin propagators of diagrams [2](b) and (c). In 
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particular, there are no contributions from mixed mesons made of one sea and one valence 
quark, so the new parameter in the mixed-action chiral Lagrangian, Amix, does not appear 
in Bk to 1-loop.^ 

We now proceed to calculate the 1-loop contributions to from Figure [2J The connected 
diagrams, [2](a), (d), and (e), combine to give the result^ 

Bo f 



+ ml.y + m\ + my 



(38) 



67r2 J (27r 

The contribution from the disconnected diagrams, [2](b) and (c) , is somewhat more tedious 
to evaluate because of the double poles in the hairpin propagators: 

M,,sc = j -^Mly + [DiAq) + Dly{q) - 2Diy{q)] , (39) 
where Dl- is the taste-singlet disconnected propagator of Eq (ITS]) . Making use of the identity 



D,, + Dyy - 2D^y = {m\ - m'yf-^-^{D,y}, (40) 

we get the following result for the disconnected piece in the "1-1-1+1" partially quenched 
theory: 



where Rj ' is the residue arising from the double pole in the disconnected propagator; Rj , 
{Mj^yj}, and {/Xj'} are defined in Sec. IIIIDI 

To get the full expression for B^ at NLO, one must combine the 1-loop contributions 
with analytic terms that arise from tree-level matrix elements of higher-order operators: 

= ^0 + I ( ^-- + ^'^- + analytic terms') . (42) 
We discuss the analytic terms in the following subsection. 



^ This cancellation of chiral logarithms containing the parameter Amix is not unique to Bk ■ It also occurs 
in other mixed action chiral perturbation theory expressions when they are written in terms of 1-loop 
PGB masses and decay constants (rather than bare parameters), such as in the case of the I — 2 tttt 
scattering amplitude [41|. 

^ We note that, although all of the integrals in this section are divergent in four dimensions, one can choose 
a suitable regulator to make them finite before their evaluation; this regulator can then be removed after 
the results have been renormalized. 
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(d) (e) 

FIG. 2: Quark flow diagram contributions to Bk at 1-loop. One external meson is a and the 
other is a K^. The two boxes represent an insertion of the Bk operator. Each box "changes" the 
valence quark flavor from d <-> s. Diagrams (a)-(c) contribute to Fig. [T]||d), diagram (d) contributes 
to Fig. [U^e), and diagram (e) contributes to Fig. [T][f). 

C. Analytic contributions to Bk at NLO 

Next-to-leading order analytic contributions to come from tree-level matrix elements 
of NLO operators. In MAXPT, such terms can come from operators of the following order 
in the power-counting scheme: 

O(p^gb), (^(«VIgb), 0{a'), 0{ml), 0{pl^^m,), 0{a'm,). (43) 

There are many such operators in the mixed action chiral Lagrangian, however, since it is 
not necessary to separate them in fits to numerical lattice data, we do not enumerate them 
all here. Instead we use symmetry arguments to restrict the possible linearly independent 
terms as in Ref. j25| . 

First we discuss those analytic contributions which are easiest to determine. One can im- 
mediately rule out contributions from operators of 0{ppQ^) because such operators contain 
four derivatives, but tree-level matrix elements contain only two fields upon which they can 
act. One can also rule out contributions from 0{a^) operators because, at tree-level, they 
would produce terms without powers of masses in them, and the chiral symmetry of the 
valence sector requires that the kaon matrix element vanishes in the chiral limit. Finally, all 
contributions from 0{a^ppQ^) operators must be proportional to mly because the deriva- 
tives must act on the two external kaons, giving a factor of which becomes m'^y when the 
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kaons are on-shell. This leads to the first analytic term: cia?m?^y. 

All of the analytic terms that are not proportional to powers of are the same as 
in the continuum partially quenched theory. The mass dependence of these terms can be 
determined using CPS symmetry {42], where C and P are the usual charge conjugation and 
parity reversal symmetries, respectively. In QCD, S corresponds to the exchange of d and 
s quarks, however, in the mixed action lattice theory, we must impose a symmetry under 
the interchange of x and y valence quarks instead: x ^ y, rrix ^ rriy. There are only two 
linearly independent 0{m?^) terms allowed by this symmetry; we choose to write them in the 
forms C2m^y oc C2{mx + triyY and c^irrix — triyY. Operators of Oip^Q^niq) only contribute 
one new hnear combination of masses: m^^Tr(Msea), where Mgea is the A^'sea x Asea sea quark 
mass matrix. 

Finally, operators of 0{a^mq) can be shown to give no new independent contributions. 
There are three possibilities for the quark mass dependence: (m^ + m^), (mj. — my), and 
Tr(Msca)- The first term, a^irrix + rny), is already included in Ci, the second, a?{mx — my), 
is forbidden by CPS symmetry, while the last term, a^Tr(Msea), vanishes for the following 
reason. The factor Tr(Msea) always comes from the operator Str(SM"''+MS"'') when S = 1, so 
the only operators that can lead to contributions of the form a^Tr(Msea) are 0{a'^) operators 
multiplying Str(SM''" + MS"!"). However, by chiral symmetry, the O(a^) operators cannot 
generate tree-level contributions to Aix- Therefore, the a^Tr(Msea) terms also vanish. 

To summarize, the following analytic terms contribute to Bk at NLO: 

[cia'^mly, C2m\y, Cz{m\ - m\f , c^mlyimfj^ + m^^ + m^^)] . (44) 

Note that we have re-expressed them in terms of meson masses, rather than quark masses, 
because those are what one measures in a lattice simulation. 



D. Next-to-Leading Order Bk Results 

In this section we present results for a "1+1+1" theory in which mu 7^ m,d 7^ ms in 
the sea sector, for a "2+1" theory in which m„ = m^ 7^ in the sea sector, and for a 
"full QCD" -like expression in which we tune the valence- valence meson masses equal to the 
taste-singlet sea-sea meson masses. 
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Bk at NLO in the 1+1+1 PQ theory is 

1 



TJ \ PQ,1 + 1 + 1 

Bo 



1 + 



[Iconn + 4+^'] + Cia' + 



C2m 



xy 



(45) 



(46) 



The connected 1-loop contribution, Iconm comes from evaluating the integral in Eq. 

/eo„n = 2mlyl{mly) - l{m\){m\ + mly) - i{ml.){ml. + mly), 
while the disconnected contribution, /^j^^^^, comes from evaluating the integral in Eq. 

« + ^r'({M|^,,}; {/.?'}) I . (47) 



rl+l+l _ 2 _™2^2_^_^ 



n_tlie_ above expressions, £ and i are integrals regulated using the standard SXPT scheme 
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(27r)4 g2 + ^2 ^ 16^ 

^^^^ ^ . J-iim') 

(27r)4 (g2 + ^2)2 15^2 ^ ^ 



(48) 
(49) 



One can completely account for lattice finite volume effects by turning the above integrals 
into sums. This yields an additive correction to the chiral logarithms 43| : 



m 



(m^) =m' ( In— + 5r(^L) ) , (ml) 



4 ^Ki{\r\mL) 



mL 



\f\ 



(50) 



(in^ + l) +Sr{mL), SrimL) = 2j2M\r\mL) (51) 



where the difference between the finite and infinite volume result is given by S[^{mL), and 
Kq and Ki are modified Bessel functions of imaginary argument. 

Finally, the residues and sets of meson masses that appear in the 1+1+1 disconnected 
contribution are defined to be: 



(52) 



(53) 
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and the mass eigenstates of the taste singlet flavor neutral PGB's in the 1+1+1 case are^ 
1 



TJIt 



(54) 



The expression in the 2+1 case is somewhat simpler: 

PQ,2+1 



1 



^1 =1 + 



+C3 



xy xy 

{rrix — rriYY 



m 



I conn 

+ ^dic] + cio' + C2m. 
+ C4(2m^^ +m|^) , 



2 

xy 



(55) 



xy 



where the connected term is the same as in the 1+1+1 case and the disconnected term is 



(57) 



When the up and down quark masses are degenerate, the flavor-neutral, taste-singlet mass 
eigenstates are: 



(58) 



The disconnected term also becomes simple enough that we choose to show it explicitly: 
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2ml 



(59) 



with 
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{mly + ml){ml^ 
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+ 



{m\-m\){m\-m\) 



(60) 



^ Strictly speaking, in the case where mu 7^ "irf, the mass eigenstates of the flavor-neutral sector are not 
the same as the physical states, tt^ and 777. Since the mixing between these two states is negligible (and 
vanishes in the isospin limit), usually one does not make a distinction between the mass eigenstates and 
the physical states. 
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Iy = Ix{X ^ Y), 



(61) 



Note that all of the sea quark dependence appears in the disconnected terms, and that the 
sum of these terms vanishes for degenerate valence quark masses. It is clear that vanishes 
when vnx = fny, but it is not immediately obvious that the other terms go to zero. However, 
it can be shown that in the limit that mx rny, Ix —>■ —Iy and thus the sum in Eq. fl59|) 
vanishes as claimed. 

Multiple definitions of the "full QCD" point appear in the mixed action XPT literature. 
This is because the mixed Ginsparg- Wilson valence, staggered sea theory has no true full 
QCD point at finite lattice spacing. Thus any choice of how to define the full QCD point 
should be made for convenience. In this paper, we consider the two cases that most closely 
resemble the full unquenched theory. 

One possible way to define full QCD for the mixed action theory with 2+1 flavors is to 
set mx = m£,j = m^o and my = rrisj. On the lattice, this is nontrivial because it requires a 
tuning of the bare valence masses in order to set the valence PGB masses to be those of the 
taste-singlet sea PGB masses. This definition has the advantage, however, that the NLO 
expression for Bx looks very much like the continuum expression, except for the analytic 
term proportional to a^. In this case, the expression for at NLO reduces to 

full" 



B 



K 



1 + 



^0 / ^^T^'^flymly 



+ cio^ + C2mly + ^ ^ + C4(2m^ + m^), (63) 

where we have used the relationship m^^ = {4:m'^y — m^)/3 which holds in this limit. ^ This 
clearly approaches the standard result as a — > 0. 

A popular, alternative way to define full QCD in MAXPT is to set the valence-valence 
meson mass equal to the pseudoscalar taste sea-sea meson mass with the same quark content. 
For perfect G-W valence quarks, this matching condition implies that the renormalized 
valence quark mass equals the renormalized sea quark mass at tree level in xPT. This is the 



Note that the coefficients Ci and C4 are different than those in Eqs. (|^5)) and ([55)) because the mass-squared 
of the taste-pseudoscalar meson and of the taste-singlet meson differ by a contribution of 0{a^). 
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matching condition most often used in mixed action lattice simulations since the Goldstone 
pion mass vanishes in the chiral limit even at finite lattice spacing. Note that although it is 
straightforward to set m™^ = m-|^ in the XPT expressions, a lattice calculation using domain 
wall fermions for the valence quarks would require a non-trivial tuning, since the coefficient 
which renormalizes the bare domain wall mass is different from that which renormalizes the 
bare staggered quark mass. The tuning for this case has been done on the MILC lattices 
by the LHP Collaboration [44] , and they find that the renormalization coefficients differ by 
around 30%. Our formula for to NLO in MA^FT with this choice of tuning (to the 
taste pseudoscalar) differs from the above tuning to the taste singlet by terms of order a^. 
Because no simplification occurs compared to the most general formula, Eq. fl45l) . we do not 
present a new expression for the tuning to the taste pseudoscalar. 

Given that the mixed action theory explicitly violates unitarity at finite lattice spacing, 
there is no a priori reason for preferring one matching condition to another, and there are 
numerous choices one could make. In the continuum limit, however, all of the matching 
choices should be identical. At fixed lattice spacing, though, the two choices mentioned 
above have their advantages and disadvantages. The advantage of matching to the taste 
singlet is that the theory is described by the full continuum QCD formula (plus an analytic 
term), but the disadvantage is that the taste singlet has the largest mass of the 16 taste pions, 
and on the MILC lattices this mass is still quite large. Matching to the taste pseudoscalar, 
the lightest of the taste pions, then has the advantage of being closer to the physical pion 
mass, but it has the disadvantage of having more complicated xPT expressions. However, 
once the explicit xPT expressions exist, this is not much of a disadvantage. In fact, given 
the complete partially quenched xPT expressions it is advantageous to not use any tuning 
at all and to take advantage of additional partially quenched data points in order to best 
constrain the unknown parameters in the chiral fit. This may not be true for all quantities, 
especially those for which xPT is not a reasonable guide. In the following numerical analysis 
we do not assume any matching condition has been chosen; we analyze the results of our 
more general partially quenched formula. 
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IV. NUMERICAL ILLUSTRATION OF MIXED ACTION CONTRIBUTIONS TO 

Bk 



We plan on carrying out a mixed action numerical simulation of Bj^ using domain wall 
fermions on the publicly available MILC improved staggered ensembles; we therefore use the 
known parameters and previous measurements on these lattices in order to obtain numerical 
error estimates for Bk using our XPT results. Currently there are two lattice spacings with 
large statistics on these ensembles, the "coarse" MILC lattices with a ^ 0.125 fm and 
the "fine" lattices with a ~ 0.09 fm. In this section we examine the modifications to the 
continuum expression for Bk that appear due to finite lattice spacing effects; these include 
both taste-breaking errors from the staggered sea sector and finite volume errors. 

The NLO expression for Bk in a mixed action theory with 2 + 1 flavors of sea quarks 
is given in Eq. flS^ . Discretization errors lead to two contributions - the shift in the 
mass-squared of the taste-singlet sea-sea meson that appears in the 1-loop disconnected 
contribution and the analytic term proportional to a? . Because we cannot a priori know 
the value of the coefficients of the analytic terms, and because the 0{a?) analytic term 
does not arise purely from taste-violating operators, we will neglect analytic terms in this 
numerical analysis of the size of taste-breaking contributions in Bk- We choose to study 
discretization errors on the a ^ 0.125 fm "coarse" MILC lattices since taste violations will 
be more pronounced than on the a ~ 0.09 fm lattices. In particular, we use the parameters 
of the ensemble with the lightest up and down sea quark masses on the smaller volume 
{L/a = 20); this ensemble has a light quark mass of amf^^ = 0.007 and a strange quark 
mass of aml^'^ = 0.05. 

In order to estimate the size of discretization errors in Bk we calculate the percent 
difference between the 1-loop contributions to Bk with and without taste-breaking: 

E^~^°°P(m r',a^Aj) - Ej^-^° °P(mr', 0) 
5],-^°°P(mri,0) 

In this expression we have set the heavier valence bare quark mass to be equal to the sea 
strange bare quark mass so that ?7 is a function of the light valence quark mass, mj'^^, and the 
taste-singlet splitting, a^Aj. The taste singlet meson is the heaviest of all of the staggered 



' ^1-loop/ val ■ ^ ' 



sea-sea mesons, and a^A/ is approximately (450 MeV)^ on the coarse lattices [ll|. Because 



the only sea-sea mesons that contribute to the Bk at 1-loop in the mixed action theory 
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FIG. 3: Percent difference between the 1-loop contributions to Bk with and without taste-breaking 
discretization errors, Eq. (j64p . as a function of valence hght quark mass. The masses and taste 
spHttings are those of the MILC coarse ensemble (a = 0.125 fm) with ami = 0.007 and anis = 0.05. 
The vertical line shows the physical light quark mass. 



are taste-singlets, this large splitting makes the effective sea quark mass considerably larger 
than a nominal light sea quark mass of rris/ 10 or m^/? would suggest. On the fine lattices 
this splitting is less, close to (280 MeV)^, which scales appropriately according to a^a^ 



and this shows that it is necessary to approach the continuum limit in order to approach the 
chiral limit in the sea sector. Note, however, that the sea quark dependence is predicted to 
be small in the non-analytic contribution to our formulas. The sea quarks only contribute 
to the disconnected hairpin diagrams in B^, and this is only around 15% of the connected 
piece at the physical point. It will, however, be necessary to study the numerical data in 
order to determine the size of the analytic contribution, as well as to test the validity of our 
XPT formulas at the physical strange quark mass. 

We plot the percent difference, Eq. as a function of valence quark mass in Fig. [3l 

setting a^Aj to be the value measured in MILC simulations on the coarse lattices [ll|. In 
this plot, the vertical line shows the location of the physical value of the average up/down 
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quark mass, mf^^'^ ~ ms/27. One can see that for larger valence light masses, r] rapidly 
goes to zero. This is to be expected as the difference between the mass-squared of purely 
valence and purely sea mesons will ultimately be negligible for sufficiently large quark masses. 
At quark masses below m™' ^ 0.01 GeV, t] begins to increase as the valence light mass 
decreases, such that as m™' — 0, the percent difference blows up rapidly. Note, however, 
that this dramatic increase does not happen until below the physical mass, so in the region 
of interest (m[^' > 0.002 GeV) the error coming from taste violations is never higher than 
0.5%. Note also that this estimate depends upon the value of the cutoff, A^, used (Fig. [3] 
uses A^ = IGeV), although any cutoff dependence can be absorbed into the analytic terms 
which are not included in this numerical analysis. Varying A^ within the range of 0.5 GeV 
to 1.5 GeV changes this picture in the relevant light quark mass range only by a nearly 
constant vertical shift at the half a percent level. 

As for the analytic terms, if all terms up to NLO were included in Eq. flMl) this ratio would 
be explicitly independent of those analytic terms present in the continuum. The remaining 
term which is proportional to has been set to zero since we do not know a priori its value. 
In our analysis, however, by choosing a non-zero value for ci, our plot in Fig. [3] would just 
shift vertically by a nearly constant amount as a function of the quark mass in the region 
of interest. If the scaling dependence of quenched domain wall fermions with various gauge 



451]), this term will give a contribution that is on 



actions is any guide (see Fig. 7 of Ref. 
the order of a few percent. 

We now repeat the above analysis, but include errors due to the finite spatial extent of the 
lattice. Such finite volume effects can be quite noticeable at the lightest sea quark masses 
available on the MILC configurations. One might imagine that the finite volume effects in 
the mixed case would not be very different than the continuum case, since the taste violating 
effects only enter through the taste singlet meson, which has a large mass. In actuality, this is 
precisely where there could be a problem, since the heavy singlet mass appears only in the sea 
sector. Partially quenched pathologies begin to appear when m™' < nf^^.^ Consequently, if 
the sea mesons are significantly heavier than the valence mesons (as they are in the mixed 
action theory) these pathologies may become more pronounced. We will see that this is, in 
fact, the case. 



See, for example, Ref. [46 1. 
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FIG. 4: Percent difference in the 1-loop contributions to Bk in finite and infinite volume, Eq. (I65p . 
as a function of valence light quark mass. The dashed curve corresponds to r]FV with a^A/ = 
(the continuum case), and the solid curve shows rjpy with a^Aj set to its value on the coarse 
(a = 0.125 fm) MILC lattices. The sea quark masses are ami = 0.007 and anig = 0.05, and the 
spatial extent of the lattice is L = 20. 



In analogy with Eq. ( !64l) . we define rjFV to be the percent difference between the 1-loop 
contribution to Bk in the mixed theory at finite volume and Bk in the mixed theory at 
infinite volume, both including discretization errors: 



B 



l-loop,FV^^val ^2 
K 



(65) 



B],-'"°''{mr\a^Aj) 

We evaluate the above expression at a spatial lattice size of L = 20; the remaining parameters 
are the same as in the previous analysis. We then plot in Fig. H] two curves - the dashed 
curve shows the percent difference in Eq. dHSD for the continuum limit, ?7i?v'(^™\ 0), while 
the solid curve shows the same percent difference with a^Aj set to its value on the coarse 
MILC lattices. Again, the vertical line indicates the physical light quark mass. 

One can see that, for the region m™^ > 0.01 GeV, the error associated with finite volume 
effects, while not negligible, is quite small, of (9(1%) or less. Although the continuum and 
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finite lattice spacing cases are different, in eacli instance finite volume effects are small. As 
the valence light quark mass drops below the sea light quark mass, the error shoots up 
rapidly, which is expected because quenching artifacts (where finite volume effects are more 
pronounced) become more noticeable in this region. As discussed above, the mixed case 
sees these quenching artifacts at a larger valence mass, since the sea mesons are heavier. 
Although one might worry about this significant difference between the continuum and the 
mixed cases for the lighter masses (the difference is rather striking at the physical mass), 
this is not a practical problem. For example, the MILC ensemble with a light sea quark mass 
of amY^ = 0.005 has a spatial length of L = 24 as opposed to 20 for the heavier sea quark 
mass ensembles. This increase in volume for the lighter quark mass was chosen by MILC to 
reduce finite volume effects for such a light sea quark mass [2]. In our planned simulations of 
Bk, the quantity which we wish to keep large is the combination m^"'''L, which as a general 
rule-of-thumb should be 4 or more to keep the finite volume effects relatively small. Note 
though, that in Fig. HJ the renormalized valence quark mass m™^ = 0.01 GeV, and with 
L=20 this corresponds to an rriTrL ^ 3, leading to a finite volume relative error of about 
1.5%. As the valence quark mass (and thus the valence pion mass) decreases, this error goes 
up rapidly. If we want to go to lighter quark masses, we will need to use the MILC lattices 
with larger volume (L=24) so that m^^'L does not become significantly smaller than 3, and 
the finite volume corrections stay below the 2% level. The key point is that simulations are 
done to the right of this "wall" in Fig. |4]at which the error explodes, so one can correct for 
finite volume errors before performing extrapolations to the continuum and physical light 
quark mass. 



V. CONCLUSIONS 



In this work we have calculated the expression for Bk in a mixed action lattice theory 
with Ginsparg- Wilson valence quarks and staggered sea quarks to next-to-leading order in 
chiral perturbation theory. We have discussed in some detail how to extend the continuum 

Although finite vohimc xPT at can be used to correctly describe more significant finite size errors, if the 
removable 1-loop finite volume corrections to Bk are around 15%, then it is likely that the remaining 
2-loop finite volume errors will still be a few percent. Such a large systematic uncertainty from finite size 
effects is unacceptable if one is aiming for an overall error of 5% in Bk- 
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calculation to the mixed action case, and we have provided expressions for both a "1+1+1" 
partially quenched theory (m„ 7^ 7^ nis) and a "2+1" partially quenched theory {rriu = 
f^d 7^ ^s), both of which reduce to the corresponding partially quenched QCD expressions 
in the continuum limit. 

It is illustrative to compare our expression for in mixed action chiral perturbation 
theory to that for other lattice theories. Four parameters are needed to describe Bk in the 
continuum: one leading order constant, Bq, and three NLO coefficients. In the case of pure 
Ginsparg- Wilson lattice fermions, the expression for B^ contains one additional coefficient 
proportional to a^. For domain- wall lattice fermions there is an additional constant, mres, 
which comes from chiral symmetry breaking due to the finite domain wall separation. This 
term simply enters B^ as an additive shift to the quark mass and can be separately measured 
in a tree-level fit to the pion mass-squared. In the mixed action lattice theory with staggered 
sea quarks and domain wall valence quarks that we have considered here, taste-symmetry 
breaking effects produce an additive shift to the sea-sea meson mass squared. This is the 
only new term that appears in the calculation of Bk with domain-wall quarks on a staggered 
sea as compared to a pure domain- wall calculation. In contrast, the expression for Bk with 
staggered valence quarks on a staggered sea contains many new parameters, each of which 
must be determined from lattice simulations and subsequently removed in order to extract 
the value of Bk in the continuum. It is interesting to note that the expression for Bk in the 
mixed G-W, staggered lattice theory is no more complicated than that for a domain-wall 
simulation in which a different value of the domain- wall separation is used in the valence and 



sea sectors; such a "mixed" domain- wall simulation was previously proposed in Ref. [47|. In 
this case the value of mres would differ in the valence and sea sectors, and the corresponding 
expression for Bk could be gotten from our expression, Eq. fHSjl . by simply making the 
replacements mres m'^^l'^^'^'^ and a'^Aj — > m^^^. Thus the taste-singlet sea-sea meson mass 
shift can be thought of as an effective mres in the sea sector, though this effective "mres" 
scales as a^, and consequently vanishes in the continuum limit. 

Finally, we have presented a numerical analysis of the resulting expressions in which we 
have examined the size of discretization errors from taste-symmetry breaking in the sea 
sector and finite volume errors for the MILC coarse (a ~ 0.125 fm, L=20 and 24) lattice 
ensembles. We find that the non-analytic taste-breaking contributions to Bk in the mixed 
action theory are around 0.5% over the range of the extrapolation and so are quite small. The 
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finite volume effects are somewliat larger for the mixed action case than in the continuum, 
but still remain at or below the 2%- level for the values of the light quark masses used to 
generate the MILC ensembles. It will of course be necessary to study the numerical lattice 
data in order to determine the size of the analytic contribution, as well as to test the validity 
of our NLO XPT formulas at the physical strange quark mass. 

A lattice calculation of Bk using domain-wall valence quarks on top of improved stag- 
gered field configurations combines the best properties of both fermion discretizations. This 
method will be competitive with other established methods for calculating Bx, and ulti- 
mately it should give a useful constraint on the CKM matrix and phenomenology. 
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